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We investigate f-electron charge susceptibility in a two-impurity Anderson model on the 
basis of Wilson's numerical renormalization group method. The f-electron charge susceptibility 
diverges logarithmically at the critical point of this model when conduction-electron bands 
exhibit particle-hole symmetry. Although the critical point disappears without the particle-hole 
symmetry, the f-electron charge fluctuation is much more enhanced near the crossover regime 
between the Kondo-Yosida singlet and intersite spin-singlet states than that in the single- 
impurity case. This result shows that charge fluctuations are enhanced owing to the competition 
between intersite and Kondo-Yosida spin singlets. A possible scenario for the enhanced residual 
resistivity near the region where the Kondo temperature becomes comparable with the Neel 
temperatures under pressure in some heavy-fermion compounds is proposed. 

KEYWORDS: two-impurity Anderson model, numerical renormalization group, residual resistivity, quantum 
critical point 



Since the first heavy-fermion compound CeAls was dis- 
covered in 1975/ many heavy-fermion f-electron systems 
have been discovered and investigated.^ Many Ce-based 
compounds show a quantum critical point (QCP) of its 
long range ordered phase under various external condi- 
tions such as pressure and magnetic fields. Near the QCP, 
anomalous properties appear in many physical quantities 
such as the diverging specific heat coefficient, the non- 
Fermi liquid (NFL) temperature dependence of the re- 
sistivity, and the enhancement in the residual resistivity. 
Furthermore, an unconventional superconductivity ap- 
pears near the QCP in many compounds.^ Among them, 
various anomalous properties near the QCP associated 
with an antiferromagnetic (AFM) long-range order have 
attracted much attention. A qualitative explanation of 
QCPs was proposed by Doniach,"^ which is competition 
between the intersite antiferromagnetic order and local 
Kondo-Yosida singlet formation. 

"Conventional" QCPs are well described by the 
Moriya-Hertz-Millis (MHM)-type critical theory of itin- 
erant magnetism, assuming the existence of quasipar- 
ticles and Landau's Fermi liquid theory.^ This means 
that Kondo-Yosida singlet formation starts even in- 
side the AFM phase, and that f-electrons form strongly 
dressed quasiparticles with a large effective mass in many 
heavy-fermion compounds. This qualitative picture is the 
basis of the MHM theory of itinerant AFM fluctuations 
in heavy-fermion systems. 

In some cases such as CeCug-KAuj,® with x — 0.1 
and YbRh2Si2,^ the conventional MHM-type critical spin 
fluctuation theory is not directly applicable, and thus al- 
ternative schemes have been proposed, in which changes 
in the Fermi surface are emphasized. -'^'''^^ In these theo- 
ries, no Kondo-Yosida singlet is formed inside the AFM 



phase. In order to explain the anomalous properties ob- 
served in these compounds, Moriya's spin fluctuation 
theory is also extended. "'^^ In spite of these theoretical 
studies, the anomalous properties of these materials re- 
main to be elucidated. 

In this Letter, we will examine the effects of intersite 
AFM interactions on the formation of heavy quasiparti- 
cles and/or Kondo-Yosida singlets in heavy-fermion sys- 
tems and show that charge fluctuations are strongly en- 
hanced by the formation of Kondo-Yosida singlet states 
in a two- impurity Anderson model (2IAM). 

This result gives a possible explanation for the en- 
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Fig. 1. (Color online) Schematic temperature- pressure phase dia- 
gram of Ce-based heavy-fermion systems in the case that po has 
a peak (a) inside the AFM phase and (b) at the QCP. TpL is the 
so-called Fermi liquid temperature. For Yb-based compounds, 
the horizontal axis should be regarded as a "negative" pressure, 
(c) Diagrammatic representation of effective potential scattering 
(double wavy line). The arrows indicate the one-particle Green's 
function, and the wavy line represents the bare potential scat- 
tering, r is the four-point vertex function. 
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hanced residual resistivity inside the AFM phase under 
pressure in CeAl2,^^ CeCusAu,^'^ and YbNi2Ge2,"^^ as 
schematicaUy shown in Fig. 1(a), because nonmagnetic 
impurity potential is enhanced by charge fluctuations/® 
leading to the enhanced residual resistivity. Interestingly, 
the Neel temperature T/v and the Kondo temperature Tk 
become comparable at the pressure where the residual 
resistivity is enhanced. Recently, similar behavior is ob- 
served in CeRu2(Gea;Sii_a;)2 at approximately x ^ 0.8.^^ 
It is also often the case that the residual resistivity is en- 
hanced at QCP (see Fig. 1(b)). However, in these com- 
pounds, the enhanced residual resistivity is observed not 
at the QCP but inside the AFM phase. This phenomenon 
has not been well recognized or understood so far, and is 
interesting because the enhanced residual resistivity can- 
not be explained by the critical AFM fluctuations near 
the QCP.18-20 

As a preliminary step, we review how to estimate, at 
least qualitatively, the residual resistivity po in heavy- 
fermion systems following the discussions in refs. 18 and 
21; we focus on the nonmagnetic impurities, which are 
the dominant sources of the residual resistivity. 

First, note that the residual resistivity strongly de- 
pends on the effective impurity scattering strength that 
the quasiparticles are subjected to. One important point 
is that the quasiparticles mainly consist of f-clectrons in 
heavy-fermion compounds. We will focus on such a situ- 
ation hereafter. Thus, it is sufficient to take into account 
only the f-electron component of the quasiparticles. The 
nonmagnetic impurity scattering strength Uq is renor- 
malized to Wq by the electron-electron interactions, as 
shown in Fig. 1(c). The residual resistivity is then given 
as po fx A'^fI'uP, with the renormalized density of states 
(DOS) at the Fermi level Np- From a detailed analysis of 
the renormalized Uq on the basis of the Fermi liquid the- 
ory in the periodic Anderson model, ^^ it is estimated that 
Wq ~ —u^{dni/dei)/Np in heavy-fermion systems for 
q — >■ 0, where rif and Cf are the f-electron number per site 
and its energy level, respectively. Thus, in heavy-fermion 
systems, the residual resistivity is enhanced when the 
f-electron charge susceptibility Xf = —{drif/def) is en- 
hanced. 

In this Letter, by using Wilson's numerical renormal- 
ization group (NRG) method, ^^ we will investigate the 
f-electron charge susceptibility xt using 2IAM as a sim- 
plest example, which includes the competition between 
the Kondo effect and the intersite exchange interac- 
tions, and also the charge fluctuations of f-electrons. The 
model itself has been well understood, ^^ together with a 
two- impurity Kondo model (2IKM),^'*"^® in the last two 
decades. There are two stable fixed points. One is the 
Kondo- Yosida singlet fixed point, where f-electron spin 
degrees of freedom are screened by those of conduction 
electrons. The other is an intersite singlet fixed point, 
where two f-electron spins form a spin-singlet state de- 
coupled from conduction electrons. Between them, there 
is a NFL fixed point when the conduction electron bands 



exhibits particle- hole (PH) symmetry. ^^ 
The Hamiltonian for 2IAM is given by 



H ^Yl ^k4aCk^ + Yl <^tfLfa<^ 



uY-fltf-^-flif»^ 



/]Vo 



^[e*''°/Lck.+h.c.] + JSi.S2, (1) 



where f^^. is an f-electron creation operator with spin 
a =t or 4, and site index a = 1 or 2, and Sq is a spin op- 
erator of an f-electron at postion xa. cj^^ is a conduction 
electron creation operator with wave number k and spin 
a. No, Ck, U and v are the number of sites, the kinetic 
energy of conduction electrons, and the Coulomb repul- 
sion and hybridization, respectively. J is a phenomeno- 
logical intersite exchange interaction between f-electrons 
and it is assumed to be antiferromagnetic, J > 0. This 
is introduced to enlarge the parameter space of 2IAM. 
One of the physical interpretations of J is, for example, 
the superexchange interaction due to the direct hopping 
processes of f-electrons. For later purposes, we define the 
distance between the two sites as _R = |xi — X2|. We 
set U — lOD [D being half the bandwidth of conduction 
electrons) and vary ef , w, and J to examine the variations 

of Xf- 

In the NRG calculation, we transform eq. (1) into a 
form of semi-infinite one-dimensional chains, and diag- 
onalize it iteratively.^^ In the two-impurity model, we 
have two conduction electron bands, namely, even and 
odd bands reflecting the mirror symmetry of the sys- 
tem.^'* We use the logarithmic discretization parameter 
A = 3 and the low-energy 1000 states are kept at each 
NRG iteration. 

First, we discuss the case of the PH-symmetric con- 
duction electron bands. We use a constant DOS for both 
the even and odd conduction electrons in the following. 
In the flnal part of this Letter, we will examine the effects 
of PH asymmetry. 

Figure 2(a) shows the total f-electron number A^f = 
J2aaifarT.fa<7) for Cf = —0.8D and V = 0.18D as a func- 
tion of J — Jc, where (• • • ) indicates the thermal average 
and Jc is the critical value of J, Jc = 0.0018215D, for 
this parameter set. At high temperatures T, N{ gradu- 
ally increases as J increases, while the variation becomes 
steeper near J — Jc as T decreases. Corresponding to 
this, Xt has a sharp peak at J = Jc as a function of J at 
low temperatures, as shown in Figs. 2(b) and (c). 

Figure 2(c) shows that the J — Jc dependence of Xf(0) 
is logarithmic, i.e., 



Xf (0) ex - log I J - Jc\ for I J - Jc\/Jc « 1, 



(2) 



at T = 0. As for the temperature dependence of Xf , Xf (^) 
just at the critical point is shown in Fig. 3 for three values 
of Cf. Note that the charge fluctuation is enhanced as 
|ef| decreases and Xf increases with decreasing T and 
diverges at T = 0. As shown in Fig. 3, the temperature 
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Fig. 2. (Color online) (a) Total f-electron number Ni vs J — Jc for 
eight different temperatures for v/D = 0.18 with PH-symmetric 
conduction electron bands, (b) f-electron charge susceptibility 
Xf(T) vs J - Jc- (c) XfiT) for J - Jc < (dotted lines) and 
J — Jc > (full linos) on logarithmic scale. 
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Fig. 3. (Color online) Temperature dependence of the f-clcctron 
charge susceptibility Xf(^) at J = Jc(ef) for e[/D = —0.4, —0.6 
and —0.8, and v/D = 0.18 with PH-symmetric conduction elec- 
tron bands. Inset: XfC^) oil logarithmic temperature scale. 



dependence of XiC^) is a-lso logarithmic: 

Xf(T)(x-logT for J = J,. (3) 

Corresponding to the diverging f-electron charge sus- 
ceptibility at J = Jc, the imaginary part of the dynam- 
ical f-electron charge susceptibility linxf{i^) has a finite 
amplitude at the frequency u! — and at J = Jci^^ 
where a finite residual entropy log\/2 remains.^* Away 
from the critical point, Imxf(a;) is proportional to lo at 
low frequency, recovering the Fermi liquid behavior. 

The above results suggest that there is a singularity in 
the charge sector through the formation of the Kondo- 
Yosida singlet between the f- and conduction electrons 
leading to a steep change in N{. The diverging Xf is 
related to the steep change in the phase shift between 
S ^ and tt/2 in the case of the PH-symmetric conduc- 
tion electrons. In order to explain this, it is essential to 
note that the numbers of f-electrons are different between 
the two fixed points: the f-electrons are more localized in 
the inter-site spin-singlet fixed point than in the Kondo- 
Yosida singlet fixed point, and the singular \J — Jc\ de- 
pendence in Xf appears between these two stable fixed 
points, i.e., at the NFL fixed point. 

On the basis of the boundary conformal field theoret- 
ical approach in 2IKM,^^ it is expected that the spin- 
singlet Cooper channel and the charge susceptibility of 
the conduction electrons diverge at the NFL fixed point. 
It is often the case that the charge susceptibility of con- 
duction electrons diverges at NFL fixed points in various 
impurity models, such as the multichannel^^ and spin 
3/2 multipolar^^ Kondo models. 

The origin of the logarithmic divergence XfC^) oc 
— log T is related to the existence of boundary opera- 
tors with the scaling dimension A = 1/2 near the crit- 
ical point of 2IAM. Here, we assume that the operator 
contents at the critical point of the 2IAM is the same 



as that in 2IKM, which is a reasonable assumption. The 
candidate is the charge vector field operator (f>c in Table 
VI in ref. 25. 4>c can be regarded as the f-electron charge 
operator and can couple with the first descendant of the 
total charge current operator Jlj^ around the NFL fixed 
point. 

Following the discussions in the multichannel Kondo 
model, ^® we obtain the logarithmic temperature depen- 
dence of the charge susceptibility. We also obtain log- 
arithmic divergence in the specific heat coefficient. As 
discussed by Johannesson et al., in the context of the 
two-channel Anderson model, '^'^ the important point is 
that the coefficient in front of the logarithm is expected 
to be scaled by 1/D rather than by 1/7V, leading to 
a small coupling constant. Then, we can safely neglect 
it in the specific heat calculation. However, in the cal- 
culation of the charge susceptibility, we cannot neglect 
it, since there is no contribution of other operators here. 
This small coefficient is the reason why we find moderate 
divergence in the f-electron charge susceptibility in Figs. 
2 and 3. 

In order to discuss a more realistic situation, we intro- 
duce PH asymmetry to the conduction electrons below. 
The importance of the PH symmetry in 2IKM and 2IAM 
has been studied in detail. •^'^'•^^ Indeed, the model lacks 
the critical point present in the PH-symmetric case. The 
two different fixed points with phase shifts 5 — Q and 
7r/2 can be connected by a smooth crossover. Although 
it is not a transition but a crossover, the nonmonotonic 
variation in Xi is expected owing to the "hidden" NFL 
fixed point. In the numerical calculations below, we set 
ek ~ VF{k — kp), kpR = tt and D = vpkp/TT, where vp 
and kp are the Fermi velocity and momentum, respec- 
tively. For these parameters, an intersite exchange in- 
teraction is induced, which is antiferromagnetic.^^ Thus, 
we set J = 0. The density of states p± for the even (-I-) 
and odd (— ) bands of conduction electrons are given as 
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Fig. 4. (Color online) f-Electron charge susceptibility Xf(^) "^^ ^ 
for ef/D = —0.8, J = 0. (a) Single impurity model with the con- 
stant DOS, and (b) two-impurity model without PH symmetry. 
For (a), the value is multiplied by the factor 2. 



p± — A±\V ± (sinfci?)/fci?], where A± is the normahza- 
tion constant. ^^ 

Figure 4 shows the v dependence of xi for eight dif- 
ferent temperatures for (a) the single-impurity Anderson 
model with a constant DOS and (b) 2IAM. As expected, 
Xf(r) for the single-impurity model shows a monotonic 
increase as v increases in (a). This is natural since as v 
increases, the charge fluctuations of f-electrons are en- 
hanced. In the case of the 2IAM, as shown in (b), Xf C^) 
develops a peak structure at low temperatures at approx- 
imately around v ~ 0.2051? in addition to the monotonic 
increase as v increases. Around v = 0.2051?, the Kondo- 
Yosida singlet starts to form as v increases, while for 
V < 0.2051? two f-electrons form a spin singlet owing to 
the induced antiferromagnetic interaction. 

As discussed above, we have demonstrated that Xf is 
enhanced near the critical point and the crossover re- 
gion between the Kondo-Yosida singlet and the intersite 
spin-singlet states. By combining this result with those in 
refs. 18 and 21, it is expected that the residual resistivity 
is enhanced owing to the Kondo-Yosida singlet forma- 
tion. Although our analysis is based on a two-impurity 
problem, we expect that the same phenomenon exists 
in the periodic lattice models for heavy-fermion systems. 
More elaborate calculations are necessary to examine the 



present scenario for the enhanced residual resistivity in 
heavy-fermion systems. 

As for the Fermi surfaces in the Anderson lattice 
model, the occurrence of the Kondo-Yosida singlet for- 
mation strongly affects the topology of the Fermi sur- 
faces inside the AFM phase. From this point, it is ex- 
pected that the residual resistivity anomaly will have a 
close relation with the changes in the Fermi surface. In 
real materials, as shown in Fig. 1(a), the system is inside 
the AFM phase at low temperatures when the residual 
resistivity shows enhancement as a function of pressure. 
In the AFM phase, the topology of Fermi surfaces in the 
Anderson lattice model is not yet fully understood and it 
is also nontrivial to regard it to be the same as that in the 
Kondo lattice model. A recent variational Monte Carlo 
study shows the existence of two different phases in the 
AFM phase. '^^ Further theoretical studies are necessary 
to clarify the relationship between the Kondo-Yosida sin- 
glet formation and the enhanced residual resistivity. 

In summary, we have discussed the f-electron charge 
susceptibility Xf(T) near the critical point of the 
two-impurity Anderson model. We have found that 
Xi{T) diverges at the critical point of this model 
with particle-hole-symmetric conduction electrons. Even 
when the conduction electron bands are not particle- 
hole-symmetric and thus the critical point smears out, 
Xf(T) is still much more enhanced than that in the 
single-impurity case around the crossover regime between 
the Kondo-Yosida singlet and intersite spin-singlet fixed 
points. These results shed light on the effects of inter- 
site correlation on heavy-fermion formation and the en- 
hanced residual resistivity in Ce- and Yb-based heavy- 
fermion compounds. 
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